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ABSTRACT 

Qualitative arguments are presented to demonstrate that the energy density 
of magnetic fields in matter accreting onto a black hole inside the marginally 
stable orbit is automatically comparable to the rest-mass energy density of 
the accretion flow. Several consequences follow: magnetic effects must be 
dynamically significant, but cannot be so strong as to dominate; outward energy 
transport in Alfven waves may alter the effective efficiency of energy liberation; 
and vertical magnetic stresses in this region may contribute to "coronal" activity. 

Subject headings: accretion, accretion disks, black holes, MHD 
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1. Introduction 

Accretion within the last stable orbit around a black hole is a very complicated 
problem. Consequently, most work on the subject has been carried out in one or another 
simplifying limit. Some (as summarized in the text by Kato et al. 1998) have studied what 
happens when the matter behaves purely hydrodynamically (i.e., without any magnetic 
effects). Others, in order to focus attention on those properties peculiar to MHD, have taken 
the limit of "force-free" magnetic fields, a limit valid when B 2 /(87r) ^> pc 2 (e.g., Blandford 
& Znajek 1977, Okamoto 1992, Ghosh & Abramowicz 1997). Still another approach has 
been to assume an arbitrary magnetic field configuration, and compute what happens when 
matter is injected into the black hole magnetosphere from a specified location at a specified 
rate (e.g., Hirotani et al. 1992). In almost every study, attention has been restricted to 
time-steady, azimuthally-symmetric situations (including, e.g., Phinney 1983 and Punsly 
1991, in which the force-free restriction is relaxed). One of the few efforts to go beyond 
these limits is the simulation by Koide et al. (1998), which permitted non-stationarity, but 
imposed azimuthal symmetry, and used arbitrary initial conditions for both the magnetic 
field and the gas pressure distribution. 

In the portion of this literature concentrating on MHD effects, the key organizing 
question has been whether these processes can tap the energy stored in black hole rotation 
(Blandford & Znajek 1977, Phinney 1983, Takahashi et al. 1990, Punsly 1998, Livio et al. 
1998, Meier 1998). 

Unfortunately, these simplifying assumptions, attractive as they are, may be artificially 
limiting. Although one can formally guarantee the validity of the force-free approximation 
by setting the accretion rate to zero (as is sometimes done), real accretion flows with 
non-zero mass accretion rates may introduce enough inertia to invalidate that assumption. 
Similarly, if continuing accretion continually brings new plasma toward the black hole, the 
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field configuration cannot be time-steady. Furthermore, simulations of MHD turbulence in 
the non-relativistic portions of accretion disks (Stone et al. 1996, Brandenburg et al. 1996) 
strongly indicate that the field is far from azimuthally-symmetric. 

It is the burden of this Letter to argue that the limitations imposed by these simplifying 
assumptions have prevented us from seeing a number of important elements in the dynamics 
of accretion onto black holes. Moreover, there has been almost no attention to the 
simple question of when these different approximations may apply. In order to get past 
these limitations, we adopt an attitude that is entirely complementary to those adopted 
previously. In order to make the smallest possible number of limiting assumptions, we 
will pay the price of relying primarily on simple scaling arguments, bolstered by one very 
schematic calculation. However, by making the minimum number of assumptions, we will 
be able to approach the question of whether the approximations that are commonly made 
are likely to be valid. The question to which we will give the greatest attention is the 
range of validity for either the purely hydrodynamic picture, or the assumption of force-free 
magnetic field structure. In addition, rather than focussing on what happens to the spin of 
the black hole, we will instead focus on the fraction of the rest-mass energy of accreting 
matter that can be released for possible radiation. 

2. The Magnetic Field Inside the Marginally Stable Orbit 

We will estimate the state of the magnetic field inside r ms , the radius of the marginally 
stable orbit, by supposing that it is "frozen" into the accreting plasma. Its value in the 
plunging region can then be scaled from its value in the disk proper by studying the flow 
lines. This view implicitly assumes that the accretion flow carries negligible net flux, that 
is, that there are hardly any field lines extending to infinity. Otherwise, the field in this 
region would be mostly due to the net magnetic flux accumulated over the black hole's 
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accretion lifetime (as argued, e.g., in Thorne et al. 1986). 

Given the matter four-velocity u a , the field evolution is conveniently described 
(Lichnerowicz 1967) by two four-vectors: E a = u^F a ^ and B a = up(*F) al3 , where F a/3 is 
the Maxwell field tensor and *F is its dual. Note that E a = (0, E) and B a = (0, B) in the 
fluid rest frame. Two of Maxwell's Equations [V • B — and Vx£ = —(l/c)dB/dt\ may 
be written as (*F)°f = 0, where we follow the usual convention and denote a covariant 
derivative by a semicolon. The definitions of E a and B a may be inverted to write *F in 
terms of the field four- vectors, so that this pair of Maxwell's equations becomes: 



U * B P _ U P B * + f^ U[x E v 



= 0, 

;/3 



where e is the completely anti-symmetric Levi-Civita symbol. Even in these circumstances, 
the plasma conductivity should easily be high enough to make E a = everywhere, i.e., to 
ensure flux-freezing. Equation 1 then becomes 

v! i B a . (i = u%B^ - u^B a - u a u^u /3 , fl B f3 . (2) 

Here the term oc B.q has been reduced by making use of the definition of B 13 and the 
expression for *F implicitly given in equation 1. From equation 2, we see that the derivative 
of B a with respect to proper time comprises three parts: response to shear (the term 
oc u?g); response to fluid density changes (the term oc u.g); and response to departures from 
free-fall (the term oc u^up-^). 

If only gravitational forces acted on matter with the angular momentum and energy 
appropriate to a circular orbit at r ms , it would accelerate inward, but initially rather slowly 
because both the first and the second derivatives of the "effective potential" (in the sense 
of Shapiro & Teukolsky 1983) are zero there. However, in a real disk, other forces can also 
act. Because the density of the disk drops sharply inward, there is an inward pressure 
gradient force (Chen & Taam 1993). In addition, the fluctuating magnetic fields that 
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remove angular momentum from matter in the bulk of the disk (Balbus & Hawley 1998) 
should work in essentially the same way here. It is likely, therefore, that plasma moves 
from the nearly-flat portion of the "effective potential" into the truly plunging region as 
a result of non-gravitational forces. Because the magnetic forces in the disk fluctuate, we 
expect that the injection of matter inside r ms will be both time-variabile and irregular as a 
function of azimuthal angle. 

However, once in the plunging region proper, one might expect gravitational forces 
to again be dominant. This assumption may be tested for self-consistency by integrating 
equations 2 with u a as given by free-fall with fixed angular momentum L and energy E (as 
evaluated at infinity, e.g, Shapiro & Teukolsky 1983). Because the metric is time-steady 
and azimuthally-symmetric, we set d/dt = d/d<p = 0. For simplicity, we consider only 
motion in the equatorial plane, i.e. u e = 0, an idealization appropriate to thin disks. With 
regard to the magnetic field, this calculation may be viewed as describing a particular small 
field loop accreted in a way unaffected by any adjacent streams. Equations 2 writen in 
Boyer-Lindquist coordinates then reduce to the single equation 



dr ~ BL dr " BL ^- 



^^L = B r _ B «3 (3) 



The field in the fluid frame is given by B^ = A"e^-Bg L , where the tensor e£ gives the 
orthonormal tetrad for a locally inertial frame (Chandrasekhar 1983), and A" is the Lorentz 
transformation from that frame to the fluid frame. 

On the other hand, the rest-mass density of matter p as measured in the fluid frame 

changes according to: 

<91np <91n(g 1//2 w r ) 

dr dr 

where — g is the determinant of the metric, so that g 1//2 gives the scale factor for the 
differential volume element. The matter density can then be expected to fall dramatically 
as matter crosses r ms because the radial speed outside r ms is ~ a(h/r) 2 v or b in a disk 



-7- 

with thickness /i«r, whereas it becomes relativistic inside r ms . Here a is the usual 
Shakura-Sunyaev (1973) dimensionless stress, and t> or f> is the velocity of a circular orbit. 

We numerically integrated equations 3 and 4 along with the equations dx a /dr = u a 
(where r is the proper time). To mimic the effect of MHD fluctuations, we took an initial 
condition in which L is reduced below the amount required for a circular orbit at r ms just 
enough to give the matter an infall velocity much smaller than the orbital velocity. The 
energy was left fixed at the energy associated with a circular orbit at r ms . 

The result is that the magnetic field in the fluid frame does increase somewhat inside 
r ms , but not by very large factors. When a = 0, B\ is virtually constant until radii very 
close to the event horizon, and increases by only a factor of two even there. Even if B^ 
is zero initially, the slow growth of the dynamical instability means that the initial orbit 
departs only slightly from circular. Consequently, shear is relatively strong compared to 
the change in density, and B^ can grow to be comparable to B$ in a fraction of an orbit 
(if the initial infall velocity is <C v or b). In fact, the diminution of B$ due to falling density 
dominates its growth due to shear only if the fluid starts out with an infall speed close to 
v or b- Otherwise, in the absence of reconnection, B^ becomes greater and greater as the 
matter takes more turns around the black hole, but diminishes somewhat from its peak in 
the final approach to the event horizon. When a = 0, the peak B^ can be as much as ~ 30 
times the initial B^. Increasing a from to 0.998M (from here on, dimensional expressions 
are written in units in which G = c = 1) does not change these results qualitatively, except 
that the peak magnitude of B^ diminishes to ~ 1.5 times the initial value. 

At the same time, the dramatic increase in radial speed quickly leads to a many order 
of magnitude decrease in p. For example, so long as the initial infall speed is <C v or b(r ms ), 
when a = (so that r ms = 6), u r = —0.1 at r ~ 4.1; when a = 0.95 (r\ ns = 1.81), the same 
speed is reached at r ~ 1.47. 



These results may now be combined to estimate, in the rest frame of plunging matter, 
the ratio between the magnetic field energy density Ub and the rest-mass density p. As 
we have just shown, Ub is likely to be similar to, or slightly greater than, the field energy 
density in the disk. In terms of the magnetic contribution to the stress ocm, the magnetic 
energy density in the disk is [B\j ' {B r B^) dj&MPd, where the quantities with subscripts d are 
measured in the disk proper. On the other hand, the rest-mass energy density is 
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where c s is the effective sound speed in the disk. Because h/r ~ c s /v or b, the ratio between 
magnetic energy density and rest-mass density in the region of plunging orbits is 
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Numerical MHD simulations of non-relativistic accretion disks (e.g., Stone et al. 1996, 
Brandenburg et al. 1996) indicate that au — en, and that B\ is rather greater than 
(B r B$)d- We have just demonstrated that in the free-fall limit, B\ is at least as large as 
B 2 d . The metric determinant ratio is generally slightly less than unity; if it is reinterpreted 
as describing the volume element, the ratio [g(r)/ ' g(r ms )] 1 / 2 = (r/r ms )[h(r)/h(r ms ) when 
the infall is contained within a thin disk centered on the black hole's equatorial plane. 
We conclude, then, that, so long as magnetic torques play an important role in driving 
accretion in the disk proper, the magnetic field must become dynamically important in the 
plunging region. As a result, it also follows that the assumption of ballistic orbits is not 
self-consistent. 

3. Consequences 



The simple estimate of equation 6 leads to a number of qualitatively important 
consequences. First, the assumption of purely hydrodynamic flow is always inappropriate; 
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continuity in B combined with the dramatic fall in matter density ensures that the magnetic 
field is dynamically significant. Simple free-fall trajectories, perhaps modified slightly by 
pressure forces, are not a good description of the streamlines; in other words, the accreting 
matter gives a significant fraction of its kinetic energy to the magnetic field. On the 
other hand, the force-free approximation is hardly better except in the extreme limit of 
zero accretion and an independently magnetized black hole, for B 2 /87r can be no more 
than comparable to p. This conclusion should be very robust because it depends only on 
flux-freezing, mass conservation, and the assumption that a« ~ a. The only way it can be 
avoided is if there is tremendous flux loss from the plasma, which would, in its own way, 
also be very interesting. 

Viewed another way, the estimate of equation 6 shows that the plunging matter inside 
r ms must do substantial work on the magnetic field, transferring much of its kinetic energy 
to it. We are then faced with the question of what happens to this energy. If it is simply 
carried into the black hole along with the plasma, the observational consequences would be 
relatively slight. Two other possibilities are both plausible and more interesting. 

To the extent that magnetic forces alter the flow, momentum and energy are carried 
along field lines. If field lines connect matter just outside r ms with matter some ways 
inside, the inertia of the plasma outside r ms may cause the field lines to rotate closer to its 
angular frequency than the considerably larger angular frequency of matter with the same 
angular momentum that has fallen to smaller radius. A torque would then be exerted in 
which angular momentum is carried back to the disk proper. The associated work done 
on the disk would ultimately be dissipated, effectively increasing the radiative efficiency 
of the disk. That the field likely takes the form of stretched loops does not diminish this 
effect: the sense of orbital shear automatically gives B r B^ the same sign no matter what 
the sign of B r is. Note that insisting on azimuthal symmetry would eliminate this sort of 
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field structure. 

How much energy and angular momentum is carried outward depends critically on, 
among other things, whether the MHD waves travel rapidly enough that waves directed 
outward in the fluid frame actually move outward in the coordinate frame. It is at this 
point in the argument that the problem of where causal disconnection sets in becomes 
central. To quantify this issue, suppose that the dispersion relation for linear waves with 
wave-vector k = +kr in the fluid frame is u 2 = v 2 k 2 , so that the phase and group speeds 
coincide. For the waves to move outward in the coordinate frame, u must stay positive 
when the wave four-momentum is transformed to that frame. By this definition, if the 
matter follows ballistic trajectories, the minimum wave speed in the fluid frame v m i n ~ \u r \ 
when a = 0; when a > 0, the ratio v m i n /\u r \ increases with both increasing a and decreasing 
r, and can be as large as ~ 10 when a = 0.998. Since B\ ~ 87rp, we expect the Alfven 
speed in the fluid frame to be ~ c, so it may be that v > v m i n in much of the inflow region, 
particularly when a is relatively small. However, a proper evaluation of this mechanism's 
efficiency clearly requires a much more complete calculation, and that is far beyond the 
scope of this paper. 

The speculation that Alfven waves may carry away significant energy seems to be at 
variance with the long-held belief that the maximum energy available for radiation is the 
binding energy of matter at the marginally stable orbit. The traditional argument has been 
that so little time is required for matter to fall from r ms to the vicinity of the event horizon 
that the relatively slow processes of kinetic dissipation have no opportunity to transform 
any more of its energy into heat, and thence into photons. However, there is a loophole in 
this argument — Alfven wave radiation is a coherent process that may transfer a substantial 
amount of energy on the dynamical time. 

The assumption of time-steadiness becomes especially limiting in this regard. If one 
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requires the field structure to be time-steady, there can be no travelling waves. However, 
when the magnetic field inside r ms remains connected to the disk (as it must, at least 
initially, if it is carried in by accreting matter), it cannot be time-steady except in some 
long-time average sense. 

If significant energy can be removed from plunging matter, carried back out to the disk, 
and dissipated just outside r ms , the effective efficiency of accretion may be greater than 
that given by the binding energy at r ms . In a Schwarzschild metric, in principle one might 
imagine that the efficiency might approach unity; in a Kerr metric, it is even imaginable 
that the outward flux of angular momentum carried by the magnetic field is so strong that 
the matter between r ms and the event horizon is put on negative energy orbits. If that 
occurs, the efficiency would become nominally greater than unity because the accumulated 
spin energy of the black hole is being tapped. In other words, there is the potential here 
for a realization of the Penrose process, made feasible by the long-range action of magnetic 
fields (Phinney 1983, Okamoto 1992, Yokosawa 1993). Both upper bounds are, however, 
somewhat unlikely to be realized in practise. Efficient removal of energy would mean that, 
at any given radial coordinate, the kinetic energy of the flow would be smaller than in 
free-fall. If this occurs through a reduction in \u r \, the relative importance of magnetic 
forces would be diminished, perhaps leading to a self-limiting of this process. Similarly, the 
"Penrose process" would likely require the orbital frequency of the matter inside r ms to be 
reduced below the orbital frequency of the matter in the disk to which it is donating its 
angular momentum. 

Whether Alfven waves can carry a significant amount of energy and angular momentum 
outwards to a "normal" accretion disk is also an issue in the theory of advection-dominated 
accretion flows (ADAFs). As described, for example, by Narayan & Yi (1995), these 
are accretion flows in which the inflow speed is close to freefall, and in which angular 
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momentum is efficiently removed to larger radii by magnetic stresses. Their mechanics 
are therefore very similar to the mechanics of accretion inside r ms as envisaged here. The 
difference here is that there is "normal" disk at radii only factors of a few outside the region 
of interest, where the energy carried outward can be dissipated. In contrast, most ADAFs 
are supposed to span a large range in radius. 

"Coronal" activity is a second possible outlet for the magnetic energy. If B 2 ~ 8irp in 
the accreting matter, but is weaker for \z\ > h, the vertical component of the VB 2 force 
would be ~ r/h times the vertical gravity of the black hole. Strong upwelling of matter 
must then result, leading directly to the creation of a significant vertical component in the 
magnetic field. Although these vertical motions are, like all the motions in this region, 
relativistic, flux loss into the corona cannot reduce the field strength by more than factors 
of a few because the inward plasma flow amplifying the field is equally rapid. 

As magnetic loops rise vertically, the shearing of their foot points can be expected to 
lead to reconnection in exactly the same fashion as expected in the non-relativistic portions 
of accretion disks (e.g., Galeev et al. 1979, Romanova et al. 1998). So much energy is 
available that the associated dissipation could be a major contributor to the "coronal" 
activity (i.e., strong hard X-ray emission) seen so often in accreting black hole systems. 
Because Compton scattering is such an efficient energy-loss mechanism, the dissipated heat 
could be radiated in less than a dynamical time: 

tcompt ~ (^) (£) ' (7) 

where \x e is the mean mass per electron and L/Le is the luminosity of seed photons in 
terms of the Eddington luminosity. In a normal plasma, m e /fx e ~ 10~ 3 , of course. 

A further consequence of the strong vertical forces that may be expected is the possible 
excitation of jets. This mechanism appears to provide an efficient way both to eject 
magnetized plasma along the rotation axis of the system, and to heat it at the same time. 
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In sum, many of the simplifying assumptions on which intuition about accretion inside 
the marginally stable orbit has been built may have misled us: Supposing the flow to 
be time-steady prevents recognition of both travelling waves and unsteady behavior like 
magnetic reconnection; supposing it to be azimuthally-symmetric forbids consideration of 
magnetic field topologies with stretched loops; supposing it to be purely hydrodynamic 
eliminates the possibility of magnetically-mediated dynamics; and supposing that the 
magnetic field is force-free rules out consideration of how infall dynamics might alter the 
field configuration. 

I am happy to acknowledge stimulating and instructive conversations with Eric Agol, 
Mitch Begelman, Roger Blandford, Jim Pringle, and Ethan Vishniac. This work was 
partially supported by NSF Grant AST-9616922 and NASA Grant NAG5-3929. 
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